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ABSTRACT 

We discuss systematically the consistency of light gluinos with data on perturbative QCD 
from deep inelastic scattering, quarkonia, jets at LEP, and the total hadronic cross-section 
in e + e~ annihilation on the Z peak and elsewhere. We demonstrate that, in addition to 
the well-known increase in the value of a s (mz) inferred from lower-energy data due to 
the slower running of ct s in the presence of light gluinos, the value of a s (mz) extracted 
from the LEP data must also be increased, as a result of including the effects of virtual 
light gluinos. The effect of these increases in other estimates of a s (mz) is to make them 
more consistent with the value extracted from the total e + e~ cross-section, which would 
otherwise appear distinctly higher. We discuss the possibility of looking for light gluinos 
in four-jet events at LEP, and their possible implications for scaling violations at HERA. 



CERN-TH.6824/93 
ACT-3/93 
CTP-TAMU-12/93 
SHEP 92/93-14 
March 1993 



1. Introduction 

It is normally assumed in discussions of supersymmetric phenomenology that all the 
sparticles have masses so large that they are beyond the physics reach of accelerators now 
running. However, there is one persistent possible exception to this general assumption, 
namely that the gluinos are very light, weighing at most a few GeV This possibility has 
been severely constrained by a number of different experimental searches, notably at the 
CERN pp collider [TJJ, in bottomonium decays, and in fixed-target experiments looking for 
metastable particles or missing energy Nevertheless, there are two possible windows 

for light gluino masses and lifetimes that do not seem to be ruled out by these analyses. 
One is for gluinos weighing between 3 and 4 GeV with lifetimes around 10 _13 s, and the 
other is for gluinos weighing 3 GeV or more and having lifetimes between about 10 -8 and 
10 ; ".s. 

Two of us (J.E. and D.V.N.) argued some time ago together with I. Antoniadis f| 
that these windows could be explored indirectly via the indirect effects of light gluinos in 
loop diagrams. Specifically, we pointed out that the consistency of low-energy determina- 
tions of a s , for example in deep inelastic scattering, quarkonium or r decays, with high- 
energy determinations of a s (mz) constrains in principle the contributions of light coloured 
particles to the renormalization group running of a s , and in particular the possible con- 
tributions of light gluinos. We showed that the LEP jet data then available disfavoured 
light gluinos, but we said that "This conclusion should be regarded as preliminary, pend- 
ing further reduction of the experimental and theoretical uncertainties". Subsequently, 
the LEP jet data and their theoretical understanding have advanced somewhat, and the 
preferred value of a s extracted from LEP jet data has increased significantly. This has 
led to renewed questioning whether light gluinos are consistent with the LEP and other 
perturbative QCD data 0,0. 

In fact, as we show later in this paper, one cannot take blindly the values of a s (mz) 
extracted from LEP jet data and compare them directly with lower-energy data extrapo- 
lated with the slower running of a s caused by the presence of light gluinos. This is because 
the LEP values of a s {mz) have been extracted from analyses using loop corrections |8| 
in which the possible effects of light gluinos have not been included. Any analysis claim- 
ing consistency of the perturbative QCD data must include light gluinos consistently in 
all loop diagrams, not just in running ot s up from the typical energies of deep inelastic 
scattering, quarkonium or r decays to the energy of LEP. 

The main purpose of this paper is precisely to furnish such a consistent analysis, 
including gluino loop effects in the analysis of LEP jet data and the total hadronic cross- 
section on the Z peak Q. We find that including the virtual effects of light gluinos increases 
the jet values of a s (mz) by about 10 percent, and also increases the value inferred from 
the total cross-section at LEP, although by less than 2 percent. The net effect of these 
increases, combined with the change in the running of a s up from lower energies, is 
that the overall consistency of perturbative determinations of a s (mz) is improved if light 
gluinos are present, since the other determinations of a s (mz) are now in better agreement 
with that extracted from the total cross-section. However, this improvement is not very 
significant, and we do not claim that the data are inconsistent with the absence of light 

f) This was not done in 0. 
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gluinos. Nevertheless, we feel that our analysis reinforces the desirability of further direct 
searches for light gluinos. One way to do this is by examining the rate and characteristics 
of four-jet events at LEP. Another is by careful analysis of scaling violations at HERA, 
where the effective value of a s would be about 10 percent higher in the presence of light 
gluinos than without them, and real gluinos might be observable at low values of the 
Bjorken scaling variable x. 



2. The Running of a s from Low Energies 

As a warm-up exercise, we first discuss the running of the low-energy a s data up 
to mz- The one-loop expression for the QCD (3 function including light gluinos is 

^9-^1 (2.1) 



4tt 



and the two-loop expression is 



where n/ is the number of flavours. 
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These expressions are to be used for momenta Q larger than the threshold for 
gluino pair production, i.e. twice the gluino mass m g i U i no , which we shall vary over the 
range 3 to 5 GeV. The most significant low-energy data are those extracted from r 
decays, charmonium spectroscopy, bottomonium decays, and deep inelastic scattering in 



fixed-target v and \i experiments pO] . The first two of these are at Q values below m g i uino , 
whilst the gluino threshold might be within the kinematic ranges of the latter two. The 
relevant average momentum transfers in the deep inelastic data are about 5 GeV for v 
scattering and 7.1 GeV for /i scattering. We list in Table 1 the values of a s {mz) extracted 



from the low-energy data, both without gluinos, as reviewed in ref. |]10| , and with gluinos 
of various different masses. We see that the effect of including gluinos is to increase the 
extrapolated value of a s {m z ) by as much as 15 percent. Since the low-energy data in 
the absence of gluinos tend to indicate a lower value of a s {mz) than that inferred from 
LEP data, this increase is promising, but we re-analyze the LEP data themselves before 
drawing any conclusions. 

f ) This is defined by 

— — —— = (3a s (Q 
dlog Q 2 
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3. 



Virtual Gluino Effects on LEP Jet Cross-sections 



In this section we shall mainly be concerned with the three-jet cross-section, which 
is given by [|ll[ 

2a s (l-y 13 ) 2 + (l-y 23 ) 2 



1 



—daf\y 13,2/23; 

C 



(3.1) 



3vr 2/132/23 

at the tree level, where Co is the leading-order contribution to the total cross-section, and 
Vij = (Pi — Vjf l m z 2 ■ The three-jet cross-section is given at the one-loop level by an 
expression of the form 
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where we have pulled out explicitly a correction term proportional to ir 2 , and the form of 
the residual correction p(?/i3, 2/23) will be discussed shortly. We assume that the infrared 
divergence in the three-jet cross-section has been regularized by combining it with the 
infrared part of the four-jet cross-section da^ /<To, which is also of order a s 2 . 

We now consider an arbitrary event shape variable X, which is given for three-jet 
final states by (2/12, 2/23), and for four-jet final states by X' 4 '(....). The corresponding 
differential cross-section is then given by 



1 da 
do dX 



dy 13 dy 23 6(l - 2/13 - VmWX - X (3) (y 13 , y 23 )) 
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and the generic correction factor, t](X), defined in ref.[|T^] is given by 

,*Xo 3) (2/13, 2/23) 
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To an extremely good approximation, as good as the accuracy of the Monte Carlo routine 
used to perform the phase space integral for the higher order corrections, the part of the 
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higher-order correction ( |3.2| ) that depends on the light fermions is given throughout the 
three-jet region by 

p Ferm \yi 3 , to) = T R Q log(y 13 y 23 ) - (3.5) 

where quarks contribute Tif/2 to Tr, and we can neglect e + e~ — > qqq'q', so that da^ = 
throughout the three-jet kinematic range. For the case of the Fox- Wolfram variable, C, the 
accuracy of this approximation can be seen by a direct comparison with the fermionic part 
of the higher-order correction displayed explicitly in |§. The key point for this analysis is 
that whereas the gluon corrections to the three jet cross section are positive the fermion- 
dependent contribution (|3.5|) is negative. This means that, when T R is increased by the 
appearance of some new species of fermion such as a light gluino, the one-loop correction 
to the three-jet cross-section is decreased by virtual fermion effects. This phenomenon is 
exhibited in the figure for the thrust variable T and the Fox- Wolfram variable C. It is 
evident, therefore, that the value of a s {mz) extracted from a fit to such a distribution 
will be increased. 

Before analysing this effect numerically, we first make some comments about the 
mathematical formulae to be used in the fits. We have argued previously |12| that the 



tt 2 terms in the one-loop correction ( |3.2|) are likely to exponentiate over most of the 
three-jet kinematic range when higher orders are calculated. We are well aware that this 
exponentiation has not been proved, but we assume it here as a working hypothesis. In 
this case, the one-loop expression becomes 

-% = I dy 13 dy 23 9{l - y 13 - y 23 )5(X - X^\y l3 ,y 23 )) ^^ 



exp {-IT 



1 + 7rp(yi3,y23) 
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The recipe for including gluinos is simply to add to p(yi 3 ,y 23 ) the expression (|3.5|) for 
pFermi w j^] 1 an ex tra contribution to Tr of 3/2, so that the differential cross-section be- 
comes da w ' 9 ', given by 

= / dy 13 dy 23 9(l - y 13 - y 23 )6(X - X^y^)) ^^ 
<Jq aJi J ctq 

17o;o7r N 
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However, there is an ambiguity in this procedure: should the exponential factor also appear 
in front of the non-7r 2 correction term p(yi3, 2/23)? No-one knows, and we consider this as a 
theoretical uncertainty in the extraction of a s (mz)- Thus, we consider in our subsequent 
numerical results two possible expressions for the gluino corrections to the differential 
cross-section for a generic three-jet variable: one without exponentiation of the gluino 
correction: 

1 da w - g - 1 da a 



a dX a dX 2tt 



^^'^^logfo^s)-!) (3-8) 



(T V 3 
and one with exponentiation of the gluino correction 
1 da w - 9 - 1 da a 



+ Y~ j ^13^230(1 - 2/13 - VXi)$(X ~ X (3) (?/i3, y 2 z)) 



o"o dX o"o rf-X' 27r 

/17« s ttW(tJ 3) (2/13,1/23) A , , 5\ 

x exp I12-J ^ [*>&<3/*v*>) - 3 J (3.9) 

where the first term on the R.H.S. in each case is the differential cross section in the 
absence of light gluinos. 

In the histograms shown in the Figure which display the effects of including light 
gluinos, the thickness of the solid line is a measure of the uncertainties in the effects, given 
by the difference between the two above expressions. 

We have studied the following three-jet variables: thrust T, oblateness O, energy- 
energy correlations EEC, the asymmetry in energy-energy correlations AEEC, the Fox- 
Wolfram variable C, the heavy jet invariant mass M#, and the quantity Mu, which 
measures the difference between the light and heavy jet invariant masses. The experimen- 



tal values of these quantities have been taken from the review p0| . The corresponding 
values of a s {m z ) extracted assuming the presence and absence of gluinos are shown in 
Table 2. In each case, the ranges quoted for the estimates with gluinos accounts for the 
above-mentioned ambiguity, and the gluino mass is assumed to be much lighter than Mz- 
We see that the inclusion of gluinos increases a s {Mz) by about 10 percent. The bottom 
row of Table 2 gives the weighted average values of a s (Mz)- For the reason given in ref. 
|T2"f , namely that the non-7r 2 correction 77 (|3.4 ) is particularly large for these variables, we 



prefer to exclude oblateness and the AEEC from the weighted average, but we do include 
the effect of exponentiation of the 7r 2 part of the correction factor discussed in [T2j . This 



has the effect of reducing the average value of a s (Mz), as shown in the last line of Table 
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2. We find that in the presence of light gluinos 



a s {M z ) = 0.124 ±.007 (3.10) 

for the value of a s [Mz) extracted from weighted averages of three-jet data on the Z 
peak. The error is obtained from the errors quoted in |TUfl added in quadrature with the 



theoretical error on the effect of the gluinos. This is to be compared with the value of 
®s( m z) — 0.113±.006 in the absence of gluinos but with tt 2 exponentiation. The penulti- 
mate line of Table 2 gives the corresponding values of a s (mz) with and withour gluinos, 
but without 7r 2 exponentiation. 



4. Total e + e Cross-section Data 

These can be used to extract a value of a s (Mz) using the next-to- leading order 



formula 13 



a tot (e + e-) = o$*(e + e 



01 / rv \ ^ 
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The available data are well known to give a value of a s [Mz) that is consistently higher 
than other determinations, though only about one standard deviation higher than the 
latest three-jet measurements on the Z peak, in the absence of light gluinos. To extract 
the value of a s {Mz) in the presence of light gluinos, we can to this order simply add three 
more effective flavours and equate 

a s (M z ) + fa^V = oOMz) + ^ (a^V 



7T \ 7T j 7T \ 7T 

where a' s (Mz) is the value which would be extracted form the total cross-section if gluinos 



were included, and a s (Mz) is the value of 0.130 ± .012 quoted in |jT0| . This yields 



ol s {Mz) = 0.132 ±.012 (4.3) 

Thus we find yet another increase in the extracted value of a s (Mz), but by less than 2 
percent this time. For the total cross-section the higher-order corrections are smaller (both 
for the gluon and fermion corrections) because of a partial cancellation by the two-loop 
correction to the two-jet cross-section. In fact it can be seen from ( |3.5|) that the effect 
of adding light gluinos actually diverges logarithmically as one goes to the two-jet region 
{Un or y 2 3 — * 0), and this divergence is cancelled by an infrared-divergent correction 
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to the two-jet cross-sect ion. The net effect is a reduction in the higher-order correction 
for the total cross-section compared with the jet event shape differential cross-sections. 

The overall consistency of the available data on a s (Mz) is therefore improved 
by the inclusion of light gluinos, in that the spread of central values is reduced from 
0.112-0.130 to 0.124-0.136, but this effect is not conclusive. The most significant effect of 
virtual gluinos is not to reconcile the low-energy and LEP jet data, whose central values 
are not brought closer together (0.124-0.136 compared with 0.112-0.121 previously). We 
suspect that this is because the three-jet data effectively measure the a s at a momentum 
scale significantly below Mz- On the other hand, including light gluinos does bring other 
determinations of a s (Mz) into better agreement with the value extracted from the total 
hadronic cross-section at the Z peak, in that a maximum possible discrepancy of 0.112 
vs 0.130 is reduced to 0.124 vs 0.132. It has been a long-standing puzzle that total cross- 
section measurements consistently give larger values of a s , although the difference has 
never been more than one standard deviation or so. Thus, there has never really been a 
problem. Nevertheless the improved agreement is interesting and may become significant 
as the statistical errors on the measurement of a s (Mz) are further reduced. Even though 
there must be less exotic solutions to the potential problem of a mismatch in the values 
of a s [Mz) extracted from different data, the fact that data on the running of a s from 
low energies to Mz cannot rule out light gluinos, together with the improved consistency 
of the LEP data when gluinos are included, re-awakens our interest in possible direct 
searches. 



5. Search for Light Gluinos in Four- Jet Events at LEP 

It was suggested several years ago || to look for light gluinos in four-jet events in 
e + e~ annihilation, produced by gluon splitting after bremsstrahlung from an initial q or 
q. Neglecting the gluino mass, which should be a good approximation for events passing 
the selection criteria of the LEP experiments, the differential cross-section is proportional 
to that for qqq'q' final states in e + e~ annihilation, where q and q' are both light flavours. 
The differential cross-section for this process can be found in |J. Note that, as mentioned 
in ref . || , there is a factor of 3 relative to the cross-section calculated for light quarks in 
ref. ||: this is composed of a factor 6 for colour, and 1/2 for the Majorana nature of the 
gluinos. 

Two of the LEP collaborations have recently published |14], pl| detailed analyses of 



f) The divergence is actually regulated by the gluino mass, but for light gluinos we can 
neglect this for energies of order Mz- 
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four-jet final states in which they make a kinematic separation of the contributions due to 
double gluon bremsstrahlung, gluon splitting to gluons, and gluon splitting to fermions. 
These studies were motivated by the wish to measure independently the colour charges 
of gluons and quarks, and the results were consistent with the Standard Model. In the 
presence of light gluinos, the value of Tr = nj/2 for quarks measured in these experiments 
would be enhanced to 

n f /2 + 3/2 (5.1) 

The data published so far are also consistent with this enhanced value. 

Discriminating between the Standard Model and the presence of light gluinos will 
require further work on at least two fronts: the QCD radiative corrections to the four- 
jet cross-section should be evaluated, and the distinctive properties of gluino jets should 
be investigated. Gluinos are expected to decay into invisible, weakly-interacting neutral 
particles that carry away missing energy. Any experiment looking at four-jet events needs 
to consider carefully whether such an energy loss could bias the jet energy determination, 
possibly shifting gluino events to lower apparent energies where there is more background, 
or even pushing gluino jets below the experimental cuts. On the other hand, gluinos in 
the longer-lived part of the light gluino window, namely with a lifetime between 10~ 8 and 
10~ 10 s, might yield events with detectable separated vertices in the four-jet region. 

We therefore think that the light gluino window could be closed (or opened) by 
determined searches at LEP. 



6. Light Gluino Effects at HERA 

If gluinos are indeed light, they could have both indirect and direct effects at 
HERA. As we have discussed above, the preferred value of a s (Mz) extracted from LEP 
and lower-energy data is increased by about 10 percent, to around 0.13, if light gluinos are 
present .[] This difference should show up clearly as enhanced scaling violations at large 
momentum transfers at HERA, due to the reduction in the value of (3. In addition to this 
indirect effect, there could also be an observable direct effect via the production of real 
light gluino pairs at small values of the Bjorken scaling variable x. Their production would 
be a higher-order QCD effect whose evaluation, as well as that of the conventional low-a; 
background, goes beyond the scope of this paper. However, we recall that, as discussed 



f ) Although it lies beyond the scope of this paper, we note that such a large value of 



a s (m.z) would endanger the high degree of consistency of LEP data |L6[ with minimal 
supersymmetric GUTs. 
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in case of four-jet events at LEP, gluino pair production events could have distinctive 
missing energies and possibly separated vertices. 

We therefore think that the light gluino window could also be closed (or opened) 
by thorough searches at HERA. 
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Table 1 

Extraction of a s (M z ) from low energy data 



Process 


a s (M z ) 


Without gluinos 


m g i U i no = 3 GeV 


mgi U i no = 4 GeV 


mgluino = 5 GeV 


r(T->^+ M -)( i£ ;p) 
r(T^+^i-) (w or i d ) 


0.121 ± .005 
0.118 ± .005 


0.140 ± .006 
0.136 ± .006 


0.137 ±0.006 
0.134 ±0.006 


0.136 ± .006 
0.132 ± .006 


J l^i T decays 


0.113 ± .005 


0.129 ± .006 


0.127 ± .006 


0.126 ± .006 


D.I.S. 


0.112 ± .004 


0.129 ± .005 


0.127 ± .005 


0.125 ± .005 
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Table 2 

Extraction of a s (M z ) from jets data at LEP 



Event Shape Variable 


ot s {M z ) 


Without gluinos 


With gluinos 


C 


0.120 (K) 
U.lz4 (DelpniJ 
0.127 (Opal) 


0.129-0.134 

n io/i n 1 on 
U. 134-11.139 

0.137-0.142 


T 


0.126 (K) 
U.lzo (Delphi) 
0.118 (L3) 
0.127 (Opal) 


0.134-0.139 
0.131-0.135 
0.126-0.130 
0.135-0.140 


EEC 


0.126 (N) 
0.129 (Delphi) 
0.128 (Opal) 


0.143-0.152 
0.147-0.156 
0.152-0.162 


M H 


0.125 (Delphi) 
0.128 (Opal) 


n 1 /II A 1 A r 

U.141-U.145 

0.133-0.137 
0.136-0.141 


M D 


U. 141 (fS J 

0.122 (Delphi) 
0.118 (Opal) 


U. 10U-U. 100 

0.130-0.134 
0.126-0.130 





0.113 (N) 
0.122 (Delphi) 
0.121 (Opal) 


122-0 127 
0.132-0.137 
0.130-0.136 


AEEC 


0.110 (N) 
0.114 (Delphi) 
0.116 (Opal) 


0.114-0.115 
0.118-0.120 
0.120-0.122 


Weighted Average 


0.120±.006 


0.132±.007 


After 7T 2 
Exponentiation 


0.113 ± .006 


0.124 ± .007 
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Figure Caption 

Histograms of the C and thrust distributions for the three-jet region, showing the lowest- 
order differential cross-section (dotted line) and the higher-order differential cross section 
without gluinos (thin solid line) and with gluinos (thick solid line). The thickness of the 
solid line is a measure of the theoretical uncertainty in the gluino effect. 



13 



